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ABSTRACT 

In a hierarchical picture of galaxy formation virialization continually transforms gravitational po- 
tential energy into kinetic energies of the baryonic and dark matter. For the gaseous component 
the kinetic, turbulent energy is transformed eventually into internal thermal energy through shocks 
and viscous dissipation. Traditionally this virialization and shock heating has been assumed to occur 
instantaneously allowing an estimate of the gas temperature to be derived from the the virial temper- 
ature defined from the embedding dark matter halo velocity dispersion. As the mass grows the virial 
temperature of a halo grows. Mass accretion hence can be translated into a heating term. We derive 
this heating rate from the extended Press Schechter formalism and demonstrate its usefulness in semi- 
analytical models of galaxy formation. Our method explicitly conserves energy unlike the previous 
impulsive heating assumptions. Our formalism can trivially be applied in all current semi-analytical 
models as the heating term can be computed directly from the underlying merger trees. Our analytic 
results for the first cooling halos and the transition from cold to hot accretion are in agreement with 
numerical simulations. 

Subject headings: cosmology: theory - early universe - galaxies: formation 



1. INTRODUCTION 

In the tr aditional hierarch ical scen ario of galaxy 

formation dRees fc Ostrikerl [l97l I White fc Rees I 



19781 IWhite fc Frenk I 1 19911 : iKauffmann et al. I HflM 
Cole et al. I I2000D . it is assumed that the infalling 
gas is shock heated to the virial temperature of the 
hosting dark matter halo near the yirial radius. How- 
ever, recently, analytical argum ent dBirnboim fe Dekel I 
l2003fl . cosmol ogical simulations (iKeres et al. I 2005) and 
obseryat i ons ([Blanton et al. I 120051: IKauffmann et al. I 
I2003afbl : ICoorav fc Milosavljevicl I2005al lb1) have all 
shown evidence that in halos below a critial mass scale 
M cr ~ 10 12 M©, infalling gas was not shock heated at 
the virial radius. Thus in addition to the traditional 
"hot mode" of accretion, halos below the critical mass 
scale also have a "cold mode" of accretion. It has 
been shown that introducing such a critical mass scale 
is adva ntageous in semi-ana l ytical galaxy formation 
model ()Cattaneo et al. I l2006t ICoorav fc Milosavlievic I 
I2005afhh . 

As the mass of a halo increases due to merging and 
accretion, the temperature of the gas inside the halo will 
also change due to turbulence and shocks which con- 
stantly transform the merging and accreting kinetic en- 
ergy into thermal energy. Since the specific kinetic en- 
ergy gained by halo during merging and accretion is given 
by the change of virial temperature, the gas heating rate 
can be computed as T = krjdT V i r /dt/("/ — 1), where 7 is 
the adiabatic index which we will take to be 5/3 as is the 
case for primordial gas. Since T V i r oc M 2 / 3 , this implies 
that the virialization heating rate is determined by the 
mass accreti on history of the halo . This point was first 
proposed bv lYoshida et al.~l (|2003h in the context of first 
structure formation. 

Our basic assumption is that due to the turbulent na- 
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ture of the virialization heating process, it is a local pro- 
cess acting on most gas particles inside the halo. This is 
very different from the popular picture that gas is virial 
heated by a spherical virial shock sitting roughly at the 
virial radius. The validity of this assumption can be seen 
from Fig. 5 of Wise & Abel (2007), which shows two- 
dimensional slices of velocity divergence for four cosmo- 
logical simulations with and without cooling. Velocity 
divergence is a good indicator of virialization shock po- 
sitions. As can be seen in that figure, in the adiabatic 
cases, these shocks mainly exist at large radii where the 
gas virializes. But when radiative cooling is included, 
the turbulence becomes highly supersonic which then in- 
creases the frequency of shock fronts in the interior of the 
halo. Even gas deep inside the halo are affected by the 
cascading network of shocks in supersonic turbulence. 
This supports our basic assumption that the potential 
energy is first transformed to kinetic energy and is ap- 
proximately shared with all the fluid in a dark matter 
halo. Under this assumption, we compute the average 
heating rate during accretion and merging. So whether 
the infalling gas can be shock heated at the virial radius 
depends on the competition between virialization heating 
rate and the local cooling rate: infalling gas can be shock 
heated at the virial radius if virialization heating rate is 
larger than the local cooling rate, otherwise infalling gas 
must flow cold to the inner halo. 

In this work, we will compute this virialization heating 
rate and compare it to the gas cooling to find the evolu- 
tion of the critical halo mass separating the cold and hot 
flow region in galaxy formation. We argue that our ap- 
proach is advantageous as it formulates galaxy formation 
in an energy conserving fashion. 

2. FORMALISM 

In the spherical c ollapse model of halo formation 
(|Gunn fc Gott Ifl972l ). a halo formed at redshift z with 
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mass M has a radius R V i r that encloses a characteristic 
overdensity of A c , which in a ACDM cosmology we will 
take to be 178 when z > 1 a nd 356/(1 + z) when z < 1 
(a more precise fit is given bv lBrvan fc Norman (1998)). 
So M is related to R V i r by M = A c ^-Rl ir p m , where 
Pm{z) — (3H 2 /8nG)Q m (l + z) 3 is the mean matter den- 
sity at redshift z. For each halo, we define a virial veloc- 
ity V^ ir = and a virial temperature T vir = ^fV 2 ir 
where /i is the mean molecular weight. From these defi- 
nitions we find that 

llrn r 2/3 A l/3 l/3„2/3 

(AT \ 2 / 3 
Mo) + 

OnV /3 /M 1/3 (±_ 

0.3/ I 178 / V0.59 



K 
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Using dz/dt = -H [Q, m (l + zf + fl A (l + z) 2 } 1 / 2 in a 
ACDM cosmology and differentiating Eq. , one finds 

2 54V3 dz 
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where dM/dz is the halo mass accretion rate, which 
in simulations can be found routinely by constructing 
halo merger trees. In this work, we will use a semi- 
analytical appr oaches to compute dM/dz pr oposed by 
Ivan den Boschl (|2002f l. IWechsler et~al~l ([20021 ) have pro- 
posed a different form of fitting formula for dM/dz, but 
their formula fitted well with van den Bosch's formula 
with redefined parameters (|van den Bosch Il2002l) . Fig.Q] 
shows the mass accretion history of halos with current 
mass 10 10 ,10 n ,2 x 10 12 , 10 13 , 10 14 , 10 15 M Q using van 
den Bosch's formula. 

3. IMPLICATIONS 
3.1. Virilization heating versus cooling 

To to see whether infalling gas can be shock heated to 
the virial temperature at the virial radius, we need to 
compare the heating rate ((2|) to the gas cooling rate at 
the virial radius. 

To compute the cooling rate, when T > 10 4 K, we use 
the tabulated cooling function fo r gas in collisional ion- 
ization equilibrium computed by ISutherland fc Dopita I 
(1993); when the gas temperature is smaller than 
10 4 K, atomic line cooling is insufficient and cooling 
rate is dominated by H2 ro-vibriational transition, we 
will use the fitting form ula for H2 cooling function in 
iRipamonti fc Abell (I2005T) assuming a universal H2 frac- 
tion f H2 = IO- 3 (|Tegmark et al. I Il997t lAbel et al. I 
1998). 

FigElshows the heating rate T and cooling rate per par- 
ticle n(r V i r )A(T V i r ) at virial radius for halos with current 
mass 10 10 , 10 11 , 2 x 10 12 , 10 13 , 10 14 , 10 15 M Q . We assume 



that the gas has an isothermal density profile and we plot 
the case for both metal-free and solar-metallicity gas. As 
can be seen from Fig. [21 in the case of metal-free gas, for 
10 10 (10 11 ) M Q halo, cooling is always at least three (two) 
orders of magnitude larger than heating; for 2 x 10 12 M© 
halo, cooling always dominates over heating but they are 
comparable around z ~ 1; for 10 13 (10 14 , 10 15 ) M halo, 
heating overtakes cooling at z ~ 3.5(5,6). In summary, 
in the case of metal free gas, for haloes with current mass 
smaller than 2 x 10 12 M Q , cooling always dominate over 
heating, while for halos with current mass larger than 
2 x 10 12 M Q , heating will overtake cooling when redshift 
is smaller than some critical redshift z cr (M) which is a 
increasing function of current halo mass. In all cases, 
the heating rate decreases rapidly when z — > due to 
the rapidly decreasing accretion rate in a dark energy 
dominated Universe. 

From Fig. [2]we can see that cooling dominates heating 
around redshift 10 while the cooling rate is still domi- 
nated by atomic line cooling. So our framework predicts 
that gas falling into halo at redshift ~ 10 generally can- 
not be shock heated at the edge of the halo if their virial 
temperature are greater than 10 4 K. 

Abel et al. (2000, 2002) showed that the first gener- 
ation of stars formed inside halos of mass ~ 10 6 M at 
redshift > 20. From Fig. [TJ we can see that this corre- 
sponds roughly to halo of current mass 10 10 M . Then 
from Fig. [21 we can see that virialization heating is com- 
parable to the cooling rate of those mi nihalos at high 
redshi ft, consistent with the conclusion of lYoshida et al. 1 
(12001 . 

3.2. Critical mass scale for shock heating at small 
redshift 

As we have seen in Fig. [51 virialization heating dom- 
inates cooling at small redshift for halos above a cirit- 
cal mass. Fig. [3] shows the evolution of the critical 
halo mass M cr where the rates are equal. We have 
considered bo th an isothermal de nsity profile and an 
NFW profile (jNavarro et al. II1997D . For NFW profile, 
we took the concentration parameter c = 12, so the 
density at the virial radius is 0.17 of the mean den- 
sity of the halo. It can be seen that for zero metallic- 
ity, M cr ~ 10 115 - 10 12 M , while for solar metallicity, 
M cr ~ 10 12 ' 5 — 10 13 , with both of them changing very 
slowly at redshift range 0.5 < z < 3. This is consis- 
tent with recent cosmological simulations which found 
that below redshift 3, only gas in halos above an al- 
most redshift independent critical mass M cr ~ 10 12 M Q 
can be shock heated at the virial radius wh ile below th is 
mass there is a cold mode of gas accretion dKeres et al 
2005tlBirnboim fc D ckcl 20031; iDekel fc Birnboim 11200 
Cattaneo et al. | |2006D. Using a spherical symmetric sta- 
bilitv analysis. iBirnboim fc Dekel I (|2003f ) explained this 
to be due to the instability of virial shock in small ha- 
los. Our result implies a complementary interpretation of 
this phenomenon. Note that M cr increases rapidly when 
z < 0.5 because halo accretion rate decrease rapidly due 
to the onset of dark energy domination. 

3.3. Implications for semi-analytical galaxy formation 

models 

A central task of semi-analytical galaxy formation 
models (SAM) is computing the evolution of hot and cold 
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gas f r actions inside da r k ma t ter halos (I White fc Rees I 
19781: IWhite fc Frenk I fl99~Ll iKauffmann et al. I 1199a 



Cole et al. 1 120001: iHernauist fc Springel 1 12003) . In tra- 
ditional SAM, cooling is treated explicitly as a dynami- 
cal process while heating is treated impulsively In SAM 
cooling is the result of competition between the local 
cooling time and dynamical time. The point where cool- 
ing time equals dynamical time is defined to be the "cool- 
ing radius" and is used to compute the amount of cold 
gas and then star formation rate. For heating, traditional 
SAM just assumes that when halo merge with each other, 
all gas that is not already cooled is s hock-heated to the 
virial temperature of t he new halo (Kauff mann et al. I 
119991: ICole et al. II2000D . We will see that if we treat 
heating also as a dynamical process, this is usually not 
the case. 

Firstly, we define a heating radius r^ as n(rh)A(T V i r ) = 
T. So gas can be shocked to the new virial temper- 
ature during merger for gas lying be tween the heating 
radius and R V i r - Fo llowing SAM (jKauffmann et aLl 
ll999tlCole et al. 1 12000). consider an isothermal gas den- 
sity profile n g {r) = M g /(ATr^m p R V i r r 2 ) where M g is the 
total gas mass, fi the mean molecular weight, then rh 
becomes 

( M g A(T mr ) V /2 (3) 



\iirfj,m p R V i r T 



In traditional SAM, it is assumed that all gas between 
the cooling radius and R vir can be heated to the virial 
temperature of the new halo in mergers. Fig. [4] shows 
the evolution of the heating radius (03) and the cooling 
radius defined by lWhite fc Frenkl (|1991l ) for halos of cur- 
rent mass 10 12 , 10 13 , 10 14 , 10 15 M . From Fig. g]we can 
see that the heating radius is always at least two times 
larger than the cooling radius for both the zero and solar 
metallicity cases. Therefore, the typical mass accretion 
rate does not supply sufficient energy to heat all of the 
gas at r > r coo i. Hence traditional SAM may signifi- 
cantly overestimated the amount of gas that can be shock 
heated during galaxy mergers. This has direct implica- 
tions for using SAM to compute thermal radiation, cos- 
mic ray acceleration, AGN feedback, or whatever process 
that depends on the amount of halo hot gas. The kinetic 
energy associated to the gas accreted in code mode just 
lost as cooling radiations. The gas that has been shocked 
heated at one redshift may not be shocked heated to the 
new virial temperature at a later merger. Their tempera- 
ture will keep to be the virial temperature corresponding 
to the time they cross the heating radius. However, note 
that this is true only when that gas is still outside the 
cooling radius, otherwise it will just cool down to the 
minimum temperature allowed by the gas cooling prop- 
erties. Thus, instead of assuming all uncooled gas to have 
the same virial temperature of the new halo, a more con- 
sist way is to compute a temperature profile using the 
virialization heating rate if what is to be calculated de- 
pends on the amount of hot gas. 

Using n g (r)A(T) — n g (r)T as the net cooling rate, one 
can estimate the change in r coo /. Thus the local cooling 
time is given by 

_ 2>k B T mr n g (r) 



wot 2[ ng -(r)A(T mr )-n g (r)T]> 
where n g (r) is the isothermal gas density profile. 



(4) 



If one assumes the density profile remains to be ap- 
proximately fixed during cooling, the gas in the halo will 
have cooled at time t out to a radius r coo i determined 
by the equality of c ooling time and halo dynamical time 
(jWhite fc Rees lfl978h 



t cool \X coal (t)*) — ^ — ^dy 
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Hence the heating corrected cooling radius is given by 
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A^nm p [\k B T mr V V ir + R vir T] 
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When r = 0, this reduces to t he usual formula of cooling 
radius ( White fc Frenk Ill991[ ). as expected. 

Fig. |4] shows the evolution of cooling radius with and 
without the heating correction term It can be seen 
that the amout of cooled gas is decreased only by a small 
amout compared to the case without heating correction. 
So virialization heating, while important for determining 
whether gas can be shock heated, plays a lesser role in in- 
fluencing gas cooling inside halos. This explains why the 
amount of cold gas pre dicted by SAM and numerical sim- 
ulatio n compares well (jYoshida et al. 1120021 : iHellv et al. I 
I200! . 

From Eqs. (01 and (O, we have 



dr coo i r CO oi w g (r coo i^A(T V i r ^ 



dt 



^cool Tig (j cool) A (7 ) 



(7) 



So combining Eqs. © and ([7]), the cooling rate 
Mcooi = dM coo i/dt = <±Trp g {r coo i)r 2 cool dr coo i/dt with heat- 
ing correction is given by 



cool 



-^fgVvirT'cool Tl g(t 'cool) A (7 ''vir) 



2Rl 



Tig ir cool) A(7i!ir ) 
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where r coo i is given by Eq. ©. 

In SAM, M coo i is a key quantity to compu te the evolu- 
tion o f cold gas and thus star formation rate (|Cole et al. I 
2000). Eq. © shows explicitly how this quantity will 
be modified in our formalism and thus can be directly 
applied to SAM. Furthermore, although we assumed T 
contains only virialization heating in this work, the fi- 
nal result Eq. ([8]) is actually general which shows how 
to consistently incorporate heating effect into SAM. For 
example, it can also be applied to discuss AGN heating 
(see e.g. Croton et al. 2006). 

In summary, we have found that gas heating is de- 
termined by the competition between cooling rate and 
heating rate while gas cooling is determined by the com- 
petition between cooling time and dynamical time. Our 
analysis suggests that just like current treatment of gas 
cooling in SAM, we should also treat gas heating as a 
dynamical and local process, especially when we try to 
compute physical quantities that relies on the amount of 
hot gas inside halos, e.g. thermal radiation from halo hot 
gas (e.g. Miniati et al. 2004, Furlanetto et al. 2005). 

4. CONCLUSIONS AND DISCUSSIONS 

We have computed a virialization heating rate which is 
directly related to the halo mass accretion history using 
van den Bosch's fitting formula. 

By comparing the virialization heating rate to the cool- 
ing rate, we find that gas can be shocked heated at the 
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virial radius only for large halos at low redshift and small 
halos at very high redshift. The critical halo mass com- 
puted in our framework agrees with recent simulations 
and other analytical arguments. 

Using the virialization heating rate, we also found that 
current SAMs may have significantly overestimated the 
amount of gas that can be shock heated. Our formalism 
provides an energy conserving remedy to this problem. 
On the other hand, gas cooling is primarily determined 
by the competition between cooling time and dynamical 
time, which explains the good fit of the cold gas amount 
in the literature by comparing SAMs to numerical simu- 
lation. 

Due to the energy-conserving nature of our formal- 
ism, it is also suitable to be used to compute quan- 
tities such as the thermal radiation from halo hot 
gas. Furthermore, since cosmic ray acceleration and 
generation of galactic magnetic fiel d are directly re- 
lated t o the amount of s h ocked gas (Waxm an fc Loeb"! 
2000t iLoeb fc Waxmanl l2000t iKeshet et al. I 120031; 
Medvedev et al. I 120051 : iPavlidou k Fields I l2006h . our 
formalism can also be applied to compute these pro- 
cesses. 

Finally, we would like to indicate that an important 
issue for our calculation is the absence of scatter in van 
den Bosch's formula as it is the averaged mass accretion 
rate in the sense of both space and time. As a space av- 



erage, it implies halos corresponding to high sigma peaks 
will have larger mass accretion rate than van den Bosch's 
formula and vice versa for halos corresponding to low 
sigma peaks. As a time average, it implies that for a sin- 
gle halo, accretion rate can become larger than van den 
Bosch's formula during major merger and smaller in the 
quiescent accretion epochs. Realizing the time-average 
nature of van den Bosch's formula may be quite impor- 
tant, since the final state of a halo may be quite different 
if one computes the heating rate using the true mass ac- 
cretion rate constructed from simulation rather than van 
den Bosch's form ula. However, w e also note that from 
Figures 2 and 3 in Ivan den Bosch 1 f2002), it's not a bad 
estimate that the scatter is roughly 0.5Af (z), almost in- 
dependent of redshift. From Eq. ([2]), we expect that this 
leads to a scatter of ~ 1.5 2 ' 3 « 1.3 in the heating rate. 
Our calculation captures some of the essential physical 
effects of the virialization heating process and is easy to 
implement in SAMs that use merger trees derived from 
N-body simulations or from Monte-Carlo techniques. 
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FlG. 1. — van den Bosch's fitting of halo mass accretion history for halos with current masses 10 10 , 10 11 , 10 12 - 3 , 10 13 , 10 14 , 10 15 M Q from 
top to bottom. 




FlG. 2. — Heating and cooling rate versus redshift. Dashed lines are heating rate for halo of current mass 10 10 , 10 11 , 2x 10 12 , 10 13 , 10 14 , 10 15 
Mq . Thick solid lines are the cooling rate at the virial radius assuming zero metallicity. Thin solid lines are the cooling rate at the virial 
radius assuming solar metallicity. 
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Fig. 4. — Evolution of heating radius (solid) and cooling radius with (dashed) and without (dot-dashed) heating correction in unit of the 
virial radius for halos with mass 10 12 ,10 13 ,10 14 ,10 15 M today (be ginning from the top-left plot and mass increases clockwise). Thick 
lines are for zero metallicity while thin lines are for solar metallicity. Note that the 10 12 Mq halo plot does not have a heating curve 
because in that case cooling dominates heating at the virial radius. 



